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ABSTRACT 

We report on the discovery of WASP-12b, a new transiting extrasolar planet with Rp\ = 
1.79 ± O.OQRj and Mp\ = 1.41 ± O.lMj. The planet and host star properties were derived 
from a Monte Carlo Markov Chain analysis of the transit photometry and radial velocity data. 
Furthermore, by comparing the stellar spectrum with theoretical spectra and stellar evolution 
models, we determined that the host star is a super- solar metallicity ([M/H]= 0.3lE5:?5), late-F 
(Teff = 6300tiJ^!^ K) star which is evolving off the zero age main sequence. The planet has an 
equilibrium temperature of Tcq=2516 K caused by its very short period orbit (P = 1.09 days) 
around the hot, 12th magnitude host star. WASP-12b has the largest radius of any transiting 
planet yet detected. It is also the most heavily irradiated and the shortest period planet in the 
literature. 

Subject headings: stars: planetary systems - techniques; radial velocities - techniques: photometric 
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1. Introduction 

Transiting extra-solar planets have provided 
tremendous information about the properties of 
planets outside our Solar System. Since 2006, a 
burst of new planet discoveries have been reported. 
We are now beginning to see the variety of ex- 
oplanets which exist in the Galaxy and to clas- 
sify them based on their properties. Furthermore, 



ST5 5BG, UK 

lOSchool of Physics, University of Exeter, EX4 4QL, UK 
Department of Physics and Astronomy, The Open Uni- 
versity, Milton Keynes, MK7 6AA, UK 

^■^Harvard-Smithsonian Center for Astrophysics, 60 Gar- 
den Street, Cambridge, MA, 02138 USA 

Observatoire de Geneve, Universite de Geneve, 51 Ch. 
des Maillettes, 1290 Sauverny, Switzerland 

^''Department of Physics, University of Warwick, Coven- 
try CV4 7AL, UK 



1 



due to the increasing number of planets being dis- 
covered and due to the detailed, multi-wavelength 
follow-up of a handful of very bright transiting sys- 
tems (e.g. HD 209458, HD 189733), we are able 
to provide strong observational tests of theoreti- 
cal models of planet formation and evolution. 

Exotic planets are particularly important be- 
cause they push the boundaries of our theo- 
retical understanding. HD 209458b, for exam- 
ple, confounded t heorists with its large radius 
(iBrown et al. 2001 ) . Since its discovery, a class of 
similar planets have been found suggesting these 
highly-irradiated, low-density planets are not rare. 
Here, we report on the discovery of a new extreme 
transiting extra-solar planet with a short orbital 
period, enlarged radius, and super-solar metallic- 
ity host star. 

In this paper, we first describe all the observa- 
tions we obtained to detect and analyse the tran- 
siting star-planet system (32). We describe the 
data analysis in ij3] where we present the proper- 
ties of the planet and its host star. Finally in SJH 
we discuss the planet properties in the context of 
current theoretical understanding of planet forma- 
tion. 



ments were obtained between 2006 November and 
2007 March. The data obtained in both sea- 
sons were processed with a custom bu ilt da ta 



reduction pipeline described in iPoUacco et al 
(l2006l) . and the resulting light curves were anal- 



ysed using a modifi ed box lea st-squares a lgorithm 
(IColher Cameron et all 120061 : iKovacs et al.il2002i ) 
to search for the planetary transit signature. 

The combined SuperWASP data showed a sig- 
nificant periodic dip in brightness with a period, 
P = 1.091 days, duration, r ~ 2.7 hours, and 
depth, 5 ~ 14 mmag. The improvement in of 
the box-shaped transit model over the flat light 
curve was = 719 and the signal-to-rednoise 
(|Pont et al.ll2006l ) was SNred = 15.6. A total of 
23 partial or full transits were captured by Super- 
WASP. 

There were no obvious objects blended with 
WASP-12 in the SuperWASP aperture, and the 
detected transit event was significant, therefore, 
WASP-12 was classed as a high priority target 
needing further study. In Figure [TJ we show the 
phase-folded light curve of the SuperWASP data, 
adopting the ephemeris resulting from the box- 
least squares analysis on the combined light curve. 



2. Observations 

2MASS J063032.79-h294020.4 (hereafter WASP- 
12) is a bright F9V star. It has been identified 
in several northern sky catalogues which pro - 
vide broad band optical I Zacharia s et al.^ 2004 ) 
and i nfra-red 2MASS magnitudes (fSkrutskie et al 
20061 ) and proper motion information. Coordi- 



nates and broad band magnitudes of the star are 
given in Table [T] 

2.1. SuperWASP Photometry 

Time series photometry of WASP-12 was ob- 
tained in the 2004 and 2006 seasons with the 
SuperWASP- N camera located on La Palma, Ca- 
nary Islands (jPollacco et al.ll2006l ). In 2004, when 
the target was first observed, only 820 photo- 
metric measurements were obtained between 2004 
August and 2004 September. However, the same 
fields were observed again in 2006 after upgrades 
to the telescope mount and instrument. Dur- 
ing the 2006 season, the target was observed in 
the field-of-view of two separate cameras, and 
a total of 5573 photometric brightness measure- 



2.2. Follow-up Multi-band Photometry 

After identification as a high priority transit 
candidate, follow-up photometry of WASP-12 was 
obtained using two additional telescopes with high 
spatial resolution ( < l"/pixel) compared to Su- 
perWASP (13.7"/pixel). We obtained observa- 
tions of WASP-12 and the surrounding region dur- 
ing the predicted time of transit to confirm that 
there are no eclipsing binaries within an arcminute 
of WASP-12 that may have caused of the transit 
signal in the SuperWASP data. WASP-12 also ap- 
pears to be a single star at the resolution provided 
by these data. The closest companion is 9" from 
the target and has a magnitude of 1^ ~ 18. En- 
couraged by this, we then obtained complete B, 
I, and z-band light curves of the transit. A de- 
tailed description of the follow-up photometry is 
described below. 

2.2.1. Tenagra Telescope photometry 

WASP-12 was observed using the fully robotic, 
Tenagra II, 0.81m f/7 Ritchey-Chretien telescope 
sited in Arizona, USA. The science camera con- 
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Table 1: Stellar Parameters for WASP- 12. The 
broad band magnitudes are obtained from the NO- 
MAD 1.0 catalogue. The stellar parameters are 
derived from our spectral synthesis of observed 
spectra of WASP- 12 (see ^XT]) . 



E 

m 0.00 
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tains a Ikxlk SITe CCD with a pixel scale of 
0.87'7pixel and a field of view of 14.8' x 14.8'. 
These photometric data were obtained as part of 
an observing program sponsored by the Las Cum- 
bres Observatory Global Telescope NetworlQ. The 
-B-band transit of WASP-12 was created from 227 
observations on two consecutive nights in March. 
The /-band light curve consists of 639 flux mea- 
surements taken over five non-photometric nights 
in 2008 February and March. 

Calibration frames obtained automatically ev- 
ery twilight were used in processing the images 
(bias subtraction and flat fielding ) with the Su- 
perW ASP data reduction pipeline ( Pollacco et al.l 
I2OO6I ). Object detection and aperture photometry, 
with a 7.5 pixel radius aperture, were performed 
on all the sta rs in the frame using DAOPHOT 
(jStetsonl \l98l\ ) run under IRAlH. The differen- 
tial photometry was derived from seven, non- 
variable comparison stars within the field which 
had V< 14 magnitude and J-H colors similar to 
the target star. For each of the filters, the fluxes 
from the comparison stars were summed, and then 
converted into a magnitude which was then sub- 
tracted from the instrumental magnitude of the 
target star. The resulting differential photometry 
of the i3-band light curve has an rms of '--^ 7 mmag. 

The /-band transit was created from data ob- 
tained in non-photometric conditions over five dif- 
ferent visits to the star (some at high airmass). 
These data were essential in confirming the tran- 
sit and refining the ephemeris, but show signifi- 
cant red noise systematics. Therefore, they were 
unsuitable for defining accurate transit parame- 
ters and were excluded when modelling the light 
curve. 

2.2.2. Liverpool Telescope photometry 

A full transit was observed using RATCam 
on the robotic 2 m Liverpool Telescope (LT) 
(jSteele et al.l 120041) on La Palma as part of the 
Canarian Observatories International Time Pro- 
gramme. A total of 614 images were taken in the 



Fig. 1. — The SuperWASP discovery light curve 
of WASP-12 phase-folded with a period, P = 
1.09142 days and epoch, Tq = 2454476.2321. 



^ www.lcogt.net 

^IRAF is written and supported by the IRAF program- 
ming group at the National Optical Astronomy Observa- 
tories (NOAO) in Tucson, Arizona. NOAO is operated by 
the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the Na- 
tional Science Foundation 
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Sloan z' band on the night of 18 February 2008. 
We employed 2x2 binning (0.27"/pixel) and used 
a 10s exposure time for all observations. The 
night was photometric, and the pointing was sta- 
ble. The target star did not drift more than a few 
pixels throughout the transit which allowed for 
very accurate differential photometry. 

The images were debiased, flat-fielded, and cor- 
rected for fringing using the standard RATCAM 
processing pipelin^. IRAF DAOPHOT was then 
used to obtain aperture photometry of the target 
and four, non- variable nearby comparison stars us- 
ing a 15 pixel radius aperture. The comparison 
stars were chosen primarily for brightness from 
within the hmited FOV (4.6' x 4.6') of the cam- 
era and are not necessarily matched in color to 
the target. The differential photometry was per- 
formed in the same way as described above for the 
Tenagra data, and the resulting out-of-eclipse light 
curve has 2.5 mmag precision. The z-band and 
-B-band transit light curves obtained via follow-up 
photometry are plotted in Figure [2] 

2.3. SOPHIE Spectroscopy 

High resolution spectroscopy of WASP- 12 was 
obtained between 2008 February 12-22 w ith the 
SOPHIE spectrograph (jBouchv et al.ll2006l ) on the 
1.93m telescope at the Observatoire de Haute 
Provence. We used the same observing pro- 
gram and instrument set-up as for other Super- 
WASP planets discovered in the northern hemi- 
sphere (IPollacco et al.ir2008HChristian et al.ll2008l : 



West et al.ll2008h . 

The weather was clear and reasonably stable 
throughout the run, although there were some 
fluctuations in seeing and transparency. There- 
fore, we adopted exposure times of 900s and 1080s 
depending on the transparency. This allowed us to 
achieve a signal-to- noise (S/N) per resolution ele- 
ment (at 5500A) of 30-45 in aU the spectra. The 
S/N values for all the spectra are listed in Ta- 
ble[2l The spectra were processed in real time with 
the SOPHIE instrument control and data reduc- 
tion software, and the RV measurements were ob- 
tained using a weighted cross-correlatio n method 
(jBaranne et al.l Il996t [Pepe et all l2005l ). To do 
this, we used a numerical mask constructed from 
the solar spectrum atlas corresponding to a G2V 



star. 

The Moon was bright over most of the run thus 
contaminating the spectra with scattered light. 
We removed its velocity s ignature (according t o 
the procedure described in IPollacco et al. ( 20081) ) 
from all the spectra. We then fit the result- 
ing cross-correlation function from each spectrum 
with a Gaussian to derive a value for the cen- 
tral RV, the fuU width half maximum (10.15-10.30 
km s^^), and the peak amplitude or Contrast 
(28.60-29.44%). The uncertainties for all RV mea- 
surements were derived using an empirical relation 
applicable to SO PHIE spectra tak e n in the high 
efficiency mode ( West et al. 2008t Bouchy et al 
in prep) . 



In this way, we obtained 21 RV measurements 
of WASP-12 over ten nights which have typical 
uncertainties of ~ 10 m s~^. The RV measure- 
ment obtained from each spectrum is listed in Ta- 
ble [5] along with the derived uncertainty. These 
data have a standard deviation of 130 m s~^, sig- 
nificantly greater than the individual uncertain- 
ties, and they vary sinusoidally when folded on 
the period derived from the transit photometry 
(Figure S]). 

Finally, any asymmetries in the line profiles 
were explored by measuring t he velocity span of 
the hne-bisector ( Grav 1988h ac c ording to the 
techique outlined in lOueloz et all mm. These 
measurements, also listed in Tabled show no cor- 
relation with radial velocity (Figure [3]), therefore 
it is unlikely that the observed RV variations are 
caused by star spots on the stellar surface or by 
blending with an eclipsing binary star in the sys- 
tem or close to the line-of-sight. 

Thus, we conclude that the observed RV varia- 
tions are caused by the gravitational influence of 
a planetary-mass object orbiting WASP-12. Fig- 
ure [4] shows a plot of the RV measurements phase- 
folded on the ephemeris given in Table [3] and over- 
plotted with the best fitting model radial velocity 
curve which is derived as described in the analysis 
section (ES]). 



^see |http: / /telescope.livjm.ac.uk/Info/Tellnst/Inst/RATCam/ 1 



2.4. Additional Spectroscopy 

The SOPHIE spectra were obtained in the high 
efficiency (HE) mode which is known to suffer from 
problems with corrections for the blaze shape. Ac- 
cording to the data products website, a residual 



Table 2: The radial velocity measurements of 
WASP-12 obtained with SOPHIE spectrograph. 
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Fig. 2. — Differential z-band (top) and Tenagra 
B-band (bottom) photometry of WASP-12 during 
the transit. An offset has been added to the B- 
band data for clarity. The data are phase-folded 
with the ephemeris given in Table [3l Overplotted 
are the best f it model transit l i ght cu rves using the 
formalism of lMandel fc Agoll (120021) applying the 
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Fig. 3. — Line-bisector velocity versus radial ve- 
locity measured from all the observed SOPHIE 
spectra. We adopt uncertainties of twice the radial 
velocity uncertainty for all bisector measurements. 
There is no correlation between these two param- 
eters indicating the radial velocity variations are 
not caused by stellar activity or line-of-sight bina- 
rity. 




Phase 

Fig. 4. — The SOPHIE radial velocity curve of 
WASP-12 phase- folded with the ephemeris given 
in Table [31 The solid line is the best model curve 
resulting from the orbital parameters of the system 
derived from the MCMC analysis in ^Xll The 
systemic RV value is shown by the dotted line. 
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blaze effect remains in the reduced spectra at the 
5% leveB These data provide exceptional radial 
velocity precision, but they are not always suitable 
for determining stellar parameters For example, 
the effective temperature is strongly constrained 
by the wings of the Ha line, and for hot stars 
with broad Ha wings, the residual instrumental 
features in the SOPHIE HE spectra lead to large 
uncertainties on this parameter. Therefore, in ad- 
dition to the 21 SOPHIE spectra which were used 
primarily to measure the radial velocity signature 
of the orbiting planet, several additional spectra of 
WASP-12 were obtained and used in deriving in- 
dependent measurements of the stellar parameters 
of the host star. The additional spectroscopic ob- 
servations are described below, and the derivation 
of the stellar parameters from all available spectra 
are described in HS.li 



2.4-1- Isaac Newton Telescope spectra 

Two individual spectra of WASP-12 were ob- 
tained on 2008 April 22 with the 2.5m Isaac New- 
ton Telescope and Intermediate Dispersion Spec- 
trograph (IDS). The longslit data were taken with 
the H1800V grating using a 1.2" slit which resulted 
in a resolution of (i? ~ 8000). A signal-to- noise of 
> 50 was achieved in both individual spectra by 
taking 900s exposures. The spectrum was centered 
at 6500A and covered the region from 6200-7000A, 
thereby providing measurements of Ha, the Li I 
doublet at 6708A and many narrow metal lines. 
Biases and lamp flats were obtained at the be- 
ginning of the night and Neon-Copper- Argon arcs 
were taken just before and just after the WASP-12 
observations. The spectra were reduced in a stan- 
dard way using the IRAF longslit package. We 
then averaged the two individual spectra and con- 
tinuum normalized the composite observation be- 
fore fitting for the stellar parameters. 

2-4-2- Telescopio Nazionale Galileo spectra 

We observed the target again on 29 April 
2008 using the high efficiency echelle spectro- 
graph, SARG, mounted on the 3.58 m Telescopio 
Nazionale Galileo (TNG) telescope. These data 
were taken as part of the Canarian Observatories 
International Time Programme- Three consecu- 
tive 1800 second exposures were taken using the 



yellow filter and grism. The spectra were binned 
2 X 1 in the spatial direction at the time of obser- 
vation to reduce the readout time. A slit width of 
0.8" was adopted which resulted in a spectral reso- 
lution of R~57000. Calibration images, including 
bias frames, lamp flat-field frames, and Thorium- 
Argon arcs, were taken at the beginning of the 
night and used in processing the target spectra 
with the REDUCE echelle dat a reduction pack- 
age ( Piskunov fc Valenti 2002f) . Special care was 
taken to provide an accurate flat-fielding of the 
data. The three individual reduced spectra were 
averaged on an order-by-order basis to produce a 
final merged spectrum which was then used in the 
determination of the stellar parameters described 
in SjSrl 

3. Analysis 

3.1. Spectroscopic analysis 

Three spectra of WASP-12 were derived from 
observations with the SARG, SOPHIE, and IDS 
spectrographs. Each independent spectrum was 
compared with synthetic spectra to determine the 
effective temperature Toff, gravity logg, metallic- 
ity, [M/H], and projected stellar rotation usini 
of WASP-12. Our spectral s ynthesis technique 
closely follows the procedure of I Valenti fc Fischer 
(|2005l ) (hereafter V F05), and a detailed d escrip- 
tion can be found in lStempels et al" (2007). 

Two additional parameters, microturbulence 
and macroturbulence, are incorporated into the 
spectral synthesis to characterise turbulent mixing 
and con vection in th e upper layers of stellar atmo- 
spheres ( Gravl [l988h . Their chosen values affect 
the derived stellar properties such that microtur- 
bulence anti-correlates strongly with metallicity, 
and macroturbulence affects the line broadening, 
and therefore the v sin i measurement. In our spec- 
tral synthesis, we closely follow VF05, so that our 
results can be compared directly with this exten- 
sive spectroscopic analysis of planet-hosting stars. 
For the microturbulence, we adopt their value of 
Vmic — 0.85 km s~^, but we note that other em- 
pirical studies of main s equence F-stars suggest 



highe r values for Vmic (see lMontalban fc D'Antona 
20071 ). For the macroturbulence, we use the em- 
pirical linear relation with temperature provided 
in VF05 to derive a value of Vmar — 4.8 km s~^. 



e | http://www.obs-hp.fr/www/guide/sophie/data_products.html | However, larger values are not excluded for hot 
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stars like WASP-12. This is due to the difficulty 
in accurately measuring macroturbulence for early 
type stars where rotational broadening typically 
dominates the line widths. Furthermore, only 79 
stars in the VF05 sample have Teff > 6200 K, so 
the empirical relation is not very well defined in 
this regime. 

We used the SARG data to derive our best 
measurement of the parameters of the host star. 
This spectrum is of high resolution and good qual- 
ity with no known residual instrumental features. 
Four spectral regions (shown in Figure [5]) were se- 
lected for the fit because they are particularly sen- 
sitive to one or more of the parameters we aim to 
derive. A simultaneous fit to these four regions of 
the spectrum yielded a Teff = 6290 K, log 5 = 4.38, 
and [M/H]=0.30. The line broadening was equiv- 
alent to that of the spectral resolution, therefore 
we were only able to derive an upper limit on the 
rotational broadening of v sini < 2.2 km s~^. This 
is derived by subtracting (in quadrature) the esti- 
mated value of macrotubulence (4.8 km s^^) from 
the width of the smallest resolvable resolution cl- 
ement (5.3 km s'^ at R=57GGG). A comparison 
of the observed and best fitting model spectrum 
is shown in Figure [51 and the final stellar parame- 
ters are listed in Table [TJ We derive uncertainties 
on these properties based on the range of values 
measured from the additional analysis of the IDS 
and SOPHIE data. 

In the SOPHIE spectrum, we simultaneously 
fit the same four spectral regions given above and 
find Toff = 6175 K, logg = 4.36, and [M/H]=0.15. 
We also measure a vsini — 4.5 km s^^. Again, 
we believe this value is only an upper limit. The 
IDS spectrum spanned the region of Ha and sur- 
rounding metal lines. In this region, there are no 
strongly gravity-sensitive features, therefore, we 
fixed the logg = 4.36, which was determined from 
the SOPHIE observations, and solved for the stel- 
lar temperature, Tcs = 6495 K, and metallicity, 
[M/H]= 0.16. The resolution was too low to mea- 
sure V sini. 

In summary, we made three independent mea- 
surements of the properties of WASP-12 by com- 
paring spectroscopic observations of the star to 
model spectra. In all three analyses, we find that 
WASP-12 is a hot, slowly rotating, metal rich, 
dwarf star. We adopt the results from the analy- 
sis of the SARG spectrum as our final values for 



the parameters of WASP-12 and the uncertainties 
on the stellar parameters from the range of values 
that were determined in the three different analy- 
ses. 

3.2. Deriving Planet and Host Star Pa- 
rameters 

The multi-band light curves and radial veloc- 
ity curve of WASP-12 were analysed simultane- 
ously in a Markov-chain Monte Carlo (MCMC) 
based routine designed specifically to solve the 
multivariate problem of transiting star-planet 
systems. The rout ine is described in detail in 



Collier Cameron et al.l (|2007( ) and iPollacco et al 
(|2008l ). The results of the box- least squares anal 



ysis to the SuperWASP photometric data (de- 
scribed in § 12. ip provide an initial estimate of 
the light curve parameters. We also initially as- 
sume an eccentricity of 0.02, a systemic RV equal 
to the mean of the velocity data, and a velocity 
amplitude derived by fitting a sinusoidal velocity 
variation to the observed RVs by minimizing x^- 
To derive a first guess for the stellar mass, we 
interpolate the super-solar metallicity (Z=0.03, 
[M/H]=0.2), zero- ag e main sequenc e temp erature- 
mass relation from Girardi et al.l (2000) at the 
stellar temperature derived in the previous sec- 
tion. Toff — 6300 K. We adopt a stellar mass of 
1.28 Mq as the initial value for this parameter. 

Via the MCMC approach, the routine repeat- 
edly adopts trial parameters until it converges on 
the set of values which produce the best model 
velocity curve and model light curves. The model 
light curves are ge nerated from the analyt ic transit 
formulas found in iMandel fc Agoj (|2002l ) (adopt- 
ing the small-planet approximation) and using 
the limb darkening coefficie nts for t he ap propri- 
ate photometric filters from IClaretl (|2000l . l2004h . 
The sum of the for input data curves with 
respect to the models is the statistic used to de- 
termine goodness-of-fit. The routine also produces 
la uncertainties on all the parameters. 

3.2.1. Evolutionary Status of the Host Star 

We ran the MCMC code initially using the Su- 
perWASP light curve, the B and z-band follow up 
photometry, and the SOPHIE radial velocity curve 
as input data sets. We ran the code without im- 
posing the main sequence prior on the overall 
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1.2 - 




Wavelength (A) 



Fig. 5. — Observed SARG spectrum (grey) overplotted by the best fitting theoretical model spectrum (solid 
black line). The top panel shows the region around the Mg b triplet (5160-5190A), the second panel is the 
region around the Na I D doublet (5850-5950 A), the third panel shows the region from 6000-62 10 A with a 
large number of metal lines, and the bottom panel is the region around the Ha line (6520-6600A). These 
regions arc modelled simultaneously with spectral synthesis to derive the parameters of the host star. The 
light grey regions of the spectrum are used to determine the continuum. Note, narrow telluric emission 
features are present in second panel at the rest wavelength of the Na I doublet feature. 
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statistic. This is not unreasonable for a late-F star 
which, according to the theoretical models, has a 
main sequence lifetime of < 1 Gyr. We then de- 
termined the evolutionary status of the host star 
using the results of this run to assess whether this 
was a reasonable assumption. 

First, we examined the lithium abundance in 
WASP-12 as a possible age indicator. In the 
SARG spectrum, there is no absorption detected 
in the Li I line located at 6708A. The IDS spec- 
trum shows a broad, very shallow absorption fea- 
ture at this position, however due to the lack of 
detection in the SARG spectrum, we suspect the 
feature is due to noise or blending with other ab- 
sorption lines (e.g. Fe I at 6707. 44A). This lack of 
Li is consistent with low levels found in old open 
clusters , like M 67 (~ 4 Gyr), for a 6300 K star 



(|Sestito fc Randichll2005h . liowever a precise age 
determination cannot be derived for the star from 
this observation. 

Next, we compare the structure and temper- 
ature of the star to the s uper-solar metallicit: 



ity 

stellar evolution models of Girardi et al.l (2000) 



to constrain the age. We use the [M/H]=0.2 
(Z=0.03) tracks which are consistent with the 
measured metallicity of WASP-12, given the un- 
certainty on this parameter. Figure [6] shows a 
modified Hcrtzprung-Russel diagram comparing 
the host star to the theoretical mass tracks and 
isochrones. Here, we plot the inverse cube root of 
the stellar density, R^/A'iy^, in solar units versus 
the stellar temperature. We choose to compare the 
data to the models in this parameter space since 

1/3 

the quantity, i?*/M*' , unlike i?, or luminosity, 
is purely observational and is measured directly 
from the light curve. In addition, it is completely 
independent of the temperature determined from 
the spectrum. 

The results of the initial MCMC run provided 
a measurement of the mean stellar density. We 

1 /3 

converted the density to R^/M^ in solar units, 
and generated the same property from the mass 
and log g values in the models. We then inter- 

1 /3 

polated in the R^/M^' — T^s plane to determine 
the mass and age for WASP-12. We interpolated 
linearly along two consecutive mass tracks to gen- 
erate an equal number of age points between the 
zero-age main sequence and the coolest point at 
the end of core hydrogen burning. We then inter- 



polated between the mass tracks along equivalent 
evolutionary points to find the mass and age from 
the models that best match the stellar properties 
derived from the MCMC code and the spectral 
snythesis. According to these particular tracks, 

1 /3 

the large value for i?*/M, indicates the star has 
evolved off the zero age main sequence, but has yet 
to reach the shell hydrogen burning stage. It is in 
a position in the diagram which give it a mass, 
AU = 1.33 ± 0.05 Mo and an age of 2.0^^ j Gyr. 
To check the accuracy and precision of this result, 
we compared the star to a second set of stellar 
evolution models. When interpolating the Z=0.03 
tracks bv iYi et al.l (j2003r ). we find a similar result. 

1/3 

The position of the star in the i?,/Af, ' —T^s gives 
an age of 2.4 Gyr and a mass of 1.38 M©. 

We investigated using the rotation period of the 
star which allows for constraining the age based on 
the expected spin-down timescale. The slow rota- 
tional period argues for an old age, however with 
only an upper limit on the v sin i we are unable 
derive an age estimate using the gyrochronology 
technique ( Earned 2007 ). 

The three age-dating techniques discussed all 
suggest WASP-12 is several Gyr old, however the 
stellar evolution models provide the only definitive 
estimate of the age. Therefore we adopt a final age 
for WASP-12 of T = 2 ± 1 Gyr. We have increased 
the uncertainty to include the error in metallicity 
and the systematics in the stellar evolution mod- 
els. 

3.2.2. Final Determination of Planet and Host 
Star Parameters 

After determining the evolutionary status of 
the host star and estimating the stellar mass from 
the evolutionary models, we ran the MCMC code 
a second time. We did not impose the main 
sequence mass-radius prior, since the star has 
evolved off the zero age main sequence. Further- 
more, we adopted an initial estimate for the stellar 
mass of 1.33 M© from the comparison to the the- 
oretical tracks. 

In addition, during the radial velocity observ- 
ing run, we attempted to detect the Rossiter- 
McLaughlin effect by taking six consecutive mea- 
surements of the target during and just after the 
expected transit. We see no evidence for this ef- 
fect in the data which confirms a small v sin i for 
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WASP-12, and is consistent with the upper hmit 
on the projected rotational velocity derived from 
the spectral synthesis. We do not currently use 
the in-transit radial velocity data to derive an in- 
dependent measurement of vsini, however we do 
not wish for this sequence of observations, com- 
prising 1/4 of our RV measurements, to bias the 
fit to the system orbital parameters. A small, but 
systematic zero-point offset in these data could af- 
fect the final results. Therefore, we treated the six 
data points taken in sequence on the same night as 
a different dataset during the fitting which allows 
them to have a slight zero-point offset compared 
to all the other data points. The resulting zero 
point offset is only 0.014 m s^^. 

The results of the final fit to the light curve 
and radial velocity curve using the SOPHIE ra- 
dial velocity data, the SuperWASP photometry, 
the Liverpool Telescope z-band, and the Tenagra 
S-band data are given in Table[3]with their la un- 
certainties. The model light curves in each band 
are overplotted on the phase-folded data in Fig- 
ure [21 and the model radial velocity curve is shown 
in Figure m In the analysis, we find that WASP- 
12b is a low density planet 0.2pj) with a mass 
approximately 40% larger than Jupiter in a 1 day 
orbit around a hot, metal rich late-F star which is 
evolving off the zero age main sequence. 

We closely examined the most fundamental fit- 
ted parameters (depth, width, and impact param- 
eter) which are used to describe the shape of the 
transit itself and find the model to be a good 
match to the data. We do note that the impact pa- 
rameter, b, which is determined by the shape of the 
ingress and egress, is difficult to measure precisely 
without exceptionally good photometry. Further- 
more, this parameter strongly influences the stel- 
lar radius and thus the planet radius through the 
depth measurement. Since the impact parameter 
has such a subtle effect on the light curve, but a 
large effect on the resulting parameters, we cal- 
culate an absolute lower limit on the planet ra- 
dius by running the MCMC code with the impact 
parameter fixed to zero (b=0 model). The b=0 
model has a greater compared to the overall 
best fitting model by 15, which is significant, but 
the effect on the light curve is very subtle. How- 
ever, when comparing the results of the b=0 model 
to the overall best fit, we find a stellar radius of 
-R*_b=o = 1-46 R0 compared to i?» = 1.57 R© and 



a planet radius of -Rpi_b=o = 1.63 Rj compared to 
i?pi = 1.79 Rj. 

4. Discussion 

WASP-12b is a unique transiting planet with 
the most apparent feature being the very large 
observed radius (1.79 Rj). The planet has a 
mean density only 24% that of Jupiter, making 
it a member of the growing class of transiting gas 
giants which all have particularly large radii. Fig- 
ure [7] shows the position of WASP-12b among the 
other published transiting planets in the mass- 
radius plane. The structure of these planets, in- 
cluding HD 209458b and WASP-lb, cannot be ex- 
plained through simple, isolated planet formation 
models, and a great deal of recent theoretical work 
has gone into determining the mechanism or mech- 
anisms causing their large sizes. It is clear that 
the external environment (stellar irradiation, tidal 
forces), the internal properties (heavy element 
abundance, clouds/hazes, day-night heat transfer, 
core mass), and the evolutio nary state (age) can 



all a. ffe ct a planet's radiu s (IGuillot fc Showman 
2OO2I : iBodenheimcr ct al.^ 120031: iFortney et al ' 



2008 



2007t Burrows et al 2 0071: iFortnev et al. ^ 
Burrows et alTf2008alibl : I Jackson et al.ll2008l) . 

WASP- 12b is the most heavily irradiated planet 
yet detected. With a host star luminosity, 
L=3.48 L0, the planet experiences an incident flux 
at its substellar point of Fp=9.03x 10^ ergs cm~^ s~ 
which is twice the stellar flux experienced by 
OGLE-TR-56b or OGLE-TR-132b, the next 



most highly irradiate d planets ([Torres et al. lEool 



Burrows et al.[ [20071 ). This intense stellar ra- 
results in an equilibrium temperature, 
2516 ((l-A)/F)"-25 K, where A is the frac- 
of absorbed flux and F is the fraction of 
the planet's surface that emits at Teq- Al though 
Jupit er has an absorbing fraction, A= 0.28 l[Tavlor 



diation 
T - 

^ eq- 

tion 



1965[ ). existing evidence suggests hot Jupiters have 
much lower albedos. For example, high precision 
optical photometry of HP 209458b gives an geo- 
metric albedo of only 4% ([Rowe et al.|[2007[ ). 

It is c lear from the detailed m o delling of hot 
Jupiters ([Bodenheimer et al. I [20031: [Fortnev et al 



20071 : [Burrows et al.|[20o'7l^ that increased incident 



stellar radiation will lead to an increase in the 
planet radius by inhibiting its contraction. This 
is a function of stellar mass and age, and less mas- 
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Table 3: WASP- 12 system parameters and la error limits derived from the MCMC analysis. 



Parameter 



Symbol 



Value 



Units 



Transit epoch (BJD) 
Orbital period 
Planet/star area ratio 
Transit duration 
Impact parameter 



To 
P 

tT 

b 



2454508. 976llEJ-5jJ]^^ 

1 noi ^94+0-000003 
h0.0002 



0.0138 
0.122 
0.36 



-0.0002 
+0.001 
-0.001 
0.05 
0.06 



days 
days 

days 
R* 



Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semimajor axis 
Orbital inclination 
Orbital eccentricity 
Longitude of periastron 



7 
a 
I 
e 

CO 



0.226troSl 
19.08451°:°^^ 

83.lt\i 
0.049l°ri^ 



km s ^ 
km s~^ 
AU 

degrees 
deg 



Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 



R* 
logfl* 
P* 



1.35: 

1.57^ 



-+0.14 
-0.14 
7+0.07 
-0.07 
4 17+0.03 
^■^ ' -0.03 

S5+° °3 



Mq 
Rq 

[cgs] 
Pe 



Planet radius 
Planet mass 

Planetary surface gravity 
Planet density 

Planet temperature {A = 0,F= 



^1) 



Rp 

Mp 

Pp 
^eq 



1.79 
1.41 
2.99 
0.24 
2516 



0.09 
-0.09 
0.10 
0.10 
-1-0.03 
-0.03 
-0.03 
-0.02 
-1-36 
36 



Rj 

Mj 

[cgs] 

PJ 
K 
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Fig. 6. — R^^lMj in solar units versus effec- 
tive temperature. WASP-12 is the large solid cir- 
cle. Evolutionary mass tracks (solid lines with 
adjacent numbers labelling the mass of that line) 
and isochrones (100 Myr (soHd), 1 Gyr (dashed), 
2.0 Gyr (do t ted), 5.0 Gyr (dot-dashed) from 
Girardi et al] (|2000t) for Z=0.03 are plotted for 
comparison to the WASP-12 parameters. 



0.1 .■■ 0.2 ,.■ 0.5 




1.0 
Mass {Mj) 



Fig. 7. — Mass versus Radius of all transiting plan- 
ets which have published masses and radii. The 
data were obtained from the exoplanet encyclope- 
dia (jhttp://exopl anet.eu/ |). WASP-12b is the solid 
black circle. Lines of constant planet density are 
overplotted (dotted) in units of Jupiter density. 



sive planets are more affected by the star's ra- 
diation in this regard. The extremely irradiated 
case of the relatively massive WASP- 12b planet 
(Mp=1.41 Mj) at an age of ~ 2 Gyr is out- 
side the range of situations p r esente d in these pa- 
pers, however. iFortnev et al.l ( 2007 ) give a radius 
of 1.248 Rj for an object similar to WASP-12b 
(1 Gyr old, 1.46 Mj), but at a distance of 0.02 AU 
from a solar luminosity star. In the context of 
the models, the incident radiation from WASP- 
12 would be equivalent to putting this simulated 
planet at a distance of ~ 0.01 AU from a solar lu- 
minosity star. Despite WASP-12b being more in- 
tensely irradiated than the simulated planet of the 
same age and mass, it is difficult to see how stel- 
lar irradiation alone could result in the observed 
radius of i?pi = 1.79 Rj (or even the absolute min- 
imum radius of i?pib=o — 1-63 Rj) in a planet as 
massive as WASP- 12b. 



Burrows et al.l (|2007r ) show how enhancing the 
atmospheric opacities of extrasolar planets also 
results in increased radii by maintaining the 
heat and entropy stored in the cores over longer 
timescales. In their model, metallicity is used as a 
proxy for atmospheric heavy element abundance 
and thus opacity. Employing an abundance of 
10 X solar allows the authors to fit the radii of 
the largest planets without adding any additional 
internal heat sources. It is very likely that WASP- 
12b has a super-solar atmospheric heavy element 
abundance, given the super-solar metallicity of the 
host star. Fu rthermore, according to the models of 



Fortnev et a l. (2008), WASP-12b should be in the 
class of highly irradiated, hot planets (pM class) 
in which high-altitude molecular hazes of TiO and 
VO (condensation temperature, T=1670 K) ab- 
sorb strongly at optical wavelengths resulting in 
larger radii. For WASP-12, measuring the pri- 
mary and secondary eclipse depths as a function 
of wavelength and obtaining precise out-of-eclipse 
photometry will allow for investigating the pres- 
ence of a high altitude absorbing population in 
the planet's atmosphere. 



Although Burrow s et al.l (|2007l ) do not calcu- 
late a model specifically for WASP-12b, they are 
able to reproduce a 1.30 Rj, 1.29 Mj, 2.5 Gyr old 
planet (OGLE-TR-56b). It will be interesting to 
see if a complete modelling of the extreme environ- 
mental conditions of WASP- 12b including stellar 
irradiation, increased heavy element abundance. 
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and high altitude hazes can produce a sufficiently 
large planet radius to match the observed value. 

If WASP-12b's large radius cannot be explained 
by increasing the atmospheric opacity or other at- 
mospheric phenomena, an additional source of in- 
ternal energy will be required to explain the ob- 
servations. Dissipation of tidal energy is a pos- 
sible contributor. The timescale for circulariza- 
tion of the WASP-12b orbit is most likely very 
short, but the best model fit to the observed ra- 
dial velocity and light curves produces a non-zero 
eccentricity of e = 0.049j;Q'o|^5. Using the formal- 
ism for th e circularization timesca le taken origi- 
n ally from Goldrei ch fc Soteil ( 1966h and provided 



Bodenheimer et al.. {200^, we find 



Tcirc = 3.2^^ Myr 

which is much shorter than the 2 Gyr age of 
WASP-12, when Qp the tidal dissipation constant 
for the planet, is given the nominal value of 10^. 

The non-zero eccentricity detection is barely a 
3cr result, however, if after further observations, 
the detection persists, either the eccentricity must 
be continuously pumped by an outer planet or 
WASP12b's tidal dissipation constant must be sig- 
nificantly larger than the typically adopted value 
(10^). Since, the presence of an additional planet 
should be detectable via long term radial velocity 
monitoring, these two scenarios are distinguish- 
able and thus this system might allow detailed 
constraints to be placed on Qp, an important, 
but difficult to measure, parameter. Additionally, 
if an second planet is causing a non-zero eccen- 
tricity, the amount of internal tidal heating that 
must be dissipated is of order 5 x 10^*^ ergs s~^. 
This is a large amount of energy (2 orders of 
magnitude greater than what is calculated by 
Bodenheimer et al.l (|2003l ) for HD 209458b) which 
would have a significant effect on the radius of the 
planet. 

WASP- 12b is the hottest and largest tran- 
siting exoplanet yet detected. It has a mass, 
Mpi=1.41 Mj, radius, i?pi=1.79 Rj, and equi- 
librium planet temperature, Tcq=2516 K. The 
planet orbits a bright late-F star with a tem- 
perature, Toff=6300 K and radius, i?*=1.57 R©. 
The host star has significantly enhanced metallic- 
ity over solar, is evolving off the zero age main 
sequence, and has an age of ^ 2 Gyr. Additional 



follow-up observations of this exciting planet will 
address some of the most pressing questions in ex- 
oplanet research, in particular what mechanism or 
mechanisms are causing the large observed radii 
of some hot Jupiters and what effect does the stel- 
lar irradition and stellar mctallicity have on the 
atmospheric and structural properties of close-in 
gas giant planets. 

The SuperWASP Consortium consists of as- 
tronomers primarily from the Queen's University 
Belfast, St Andrews, Keele, Leicester, The Open 
University, Isaac Newton Group La Palma and In- 
stituto de Astroffsica de Canarias. The Super- 
WASP Cameras were constructed and operated 
with funds made available from Consortium Uni- 
versities and the UK's Science and Technology Fa- 
cilities Council. SOPHIE observations have been 
funded by the Optical Infrared Coordination Net- 
work. The data from the Liverpool and NOT tele- 
scopes was obtained under the auspices of the In- 
ternational Time of the Canary Islands. We ex- 
tend our thanks to the staff of the ING and OHP 
for their continued support of Super WASP-N and 
SOPHIE instruments. 

The radial velocity observations have been 
funded by the Optical Infrared Coordination net- 
work (OPTICON), a major international collabo- 
ration supported by the Research Infrastructures 
Programme of the European Commissions Sixth 
Framework Programme. 

This paper was based in part on observations 
made with the Italian Telescopio Nazionale Galileo 
(TNG) operated on the island of La Palma by 
the Fundacin Galileo Galilei of the INAF (Istituto 
Nazionale di Astrofisica) at the Spanish Observa- 
torio del Roque de los Muchachos of the Instituto 
de Astrofisica de Canarias 

The Liverpool Telescope is operated on the is- 
land of La Palma by Liverpool John Moores Uni- 
versity in the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de 
Canarias with financial support from the UK Sci- 
ence and Technology Facilities Council. 
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